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ABSTRACT 



We examine the wavelength dependence of flux ratios for six gravitationally lensed 
quasars using K and V images obtained at the Gemini North 8m telescope. We select 
lenses with source redshifts 2:5 < 2.8 so that K-band images probe rest-frame optical 
emission from accretion disks, while L'-band images probe rest-frame near-infrared flux 
emitted (in part) from the more extended surrounding torus. Since the observations 
correspond to different source sizes, the K and flux ratios are sensitive to structure 
on different scales and may be useful for studying small-structure in the lens galaxies. 
Four of the six lenses show differences between K and V flux ratios. In HE 0435—1223, 
SDSS 0246-0825, and HE 2149-2745 the differences may be attributable to known 
microlensing and/or intrinsic variability. In SDSS 0806+2006 the wavelength depen- 
dence is not easily attributed to known variations, and may indicate the presence of 
substructure. By contrast, in Q0142-100 and SBS 0909+523 the K and U flux ratios 
are consistent within the uncertainties. We discuss the utility of the current data for 
studying chromatic effects related to microlensing, dust extinction, and dark matter 
substructure. 



1. Introduction 

While the cold dark matter (CDM) paradigm for structure formation successfully describes 
cosmological observations on large (CMB and cluster) scales, there is notable disagreement with 
small-scale observations. Among other issues, A^-body simulations (e.g., Via Lactea - Diemand et 
al. 2008; Aquarius - Springel et al. 2008) predict the existence of numerous CDM subhalos, with 
masses M ^ 10^-10^ M©, embedded in galaxy-scale dark matter halos. This has proved troubhng 
observationally, because there are many fewer dwarf galaxies in our own Milky Way than predicted 



by CDM. Since the discrepancy may be due to baryon stripping from subhalos (e.g., [Madau et al. 



2008; Maccio et al. 2010), we need ways to probe dark matter substructure directly, regardless of 



the presence of baryonic material. 
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Gravitational lensing provides a unique way to detect CDM substructure in distant galaxies 
(e.g., |Metcalf fc Madau||200H |Ctiiba||2002| [Dalai fc Koclianek||2002D . Stars and CDM substructure 
perturb the lens potential on micro- to milli-arcsecond scales, which can have dramatic effects on 
the properties of lensed images. Most notably, lensing from stars and dark matter substructure can 



alter the flux ratios from those of smooth mass distributions. As shown by Dobler &; Keeton (2006) 



lens flux ratios depend on the size of the source compared to the size of the perturber. When the 
source is very small, it is effectively a point source for the lens substructure, so it feels an additional 
magnification boost (Sgub = Msub/^smooth) over and above the macroscopic lens properties. As the 
source increases in size, Sg^b may increase or decrease depending on the source's location relative 
to the substructure and the parity of the "macro" magnification. As the source grows still larger, 
it ceases to be affected by the substructure and Sgub 1- This phenomenon implies that by 
measuring flux ratios at different wavelengths, corresponding to different source sizes, we may be 



able to map substructure on a variety of scales (also see Moustakas & Metcalf 2003). 



Heuristically, a quasar emitting region of size Rs is significantly affected by a subhalo with 
Einstein radius Re only if Rs < Re- For typical lens and source redshifts (e.g., zi = 0.5 and 
Zs = 2.0), the Einstein radius of a subhalo of mass M is Re ^ 10^^ cm (M/Mq)^/'^. Since the 
optically emitting regions of quasars have Rs ^ 10^^-10^^ cm (Wyithe et al. 2000; Pooley et al. 



2007; Morgan et al. 2010), optical lens flux ratios are sensitive to both microlensing by stars and 



millilensing by CDM substructure. By contrast, the more extended infrared emitting regions with 
i?^ > 1 pc dChiba et aLl|2005t [Minezaki et~aLl|2009t |Agol et al.||2009D can only be affected by 



relatively massive subhalos. Comparing lens flux ratios at different wavelengths therefore makes it 
possible to constrain the amount of micro- and milli-lensing present in the system, as well as the 
sizes of the perturbers. 

Previous studies have used mid-IR observations to probe rest-frame IR emission, yielding 
evidence for subhalos with masses > IO^Mq in the lenses B1422+231 and MG 0414+0534 (|Chiba 



et all 120051 iMinezaki et al.l 120091), constraints on the mass of a luminous companion galaxy in 
H1413+117 dMacLeod et al.|2009J p1and null detections in several other systems ( |Chiba et al.|2005| 



Minezaki et al.|[2009t [Agol et "aL||2009| ). Here we extend the study of wavelength-dependent flux 
ratios by using K (2.2 /xm) and (3.8 fim) images of six lenses obtained with Gemini North during 
the 2008B semester. For source redshifts Zg < 2.8, the L'-band images correspond to rest-frame 
emission at > 1 /im where ^20-100% of the flux is thermal emission from the inner dusty torus 
( |Rowan-Robinson||1995] [Minezaki et al.||2004t [Nenkova et"aL||2008t |Agol et al.||2009D . By contrast, 
the X-band flux comes mostly from the smaller accretion disk. Thus, comparing K and V flux 
ratios may provide a sufficient source size baseline to identify substructure. 

To be sure, there are phenomena besides CDM substructure that may cause lens flux ratios 
to vary with wavelength. Optical observations probe rest-frame UV emission from the accretion 



disk of the AGN, so they are sensitive to microlensing by stars in the lens galaxy (e.g., Chartas 



^ More et al 



(2009) also constrained the mass of a luminous satellite, but using radio observations. 
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et al. 2009; Morgan et al. 2010| ). Optical flux ratios can also be altered by differential extinction 



from dust in the lens galaxy (e.g., EKasdottir et al. 2006). Finally, intrinsic variability in the 



source coupled with the lens time delay can cause flux ratios to vary with time, and the variability 
may be chromatic]^ All three effects should be attenuated at near-IR wavelengths, because the 
effective size of the accretion disk is larger and some of the flux originates in the larger dusty torus, 
but they may not be entirely absent. In particular, the importance of microlensing and intrinsic 
variability for L'-band observations will depend on the relative strengths of the accretion disk (0.01- 
0.1 pc, 10~^-10~^ arcsec) and dusty torus (0.5-5 pc, 10~^-10~^ arcsec) emission. Any significant 
wavelength-dependence in lens flux ratios is interesting, whether related to CDM substructure or 
not, so at minimum our observations can highlight systems that warrant further study. 

Most lensing constraints on CDM substructure have relied on the identification of flux ratio 
anomalies, which is best done using four-image lenses. The reason is that identifying such anomalies 
requires either well-constrained lens models (e.g., Dalai fc Kochanek||2002 ; Metcalf fc Zhao|2002 ) or 



universal magnification relations that apply only to certain four-image configurations (Keeton et al. 



2003, 2005). However, the search for wavelength dependence in flux ratios is a purely observational 



task that may provide model-independent evidence for substructure. We are therefore able to 
analyze two-image lenses for the presence of substructure for the first time. Some care is needed 
when interpreting the results (see but the observations are still valuable. 



2. Observations 



To explore the multi-wavelength behavior of flux ratios, we obtained the first K (2.2 /im) and 
(3.8 fim) images for six gravitational lenses using the Near Infrared Imager (NIRI) located at the 
Gemini North facility. Our lens sample was selected from the set of known quasar lenses that can 
be found in the CASTLES database^ Our targets (Table [l]) were chosen to have relatively bright 
{rriH < 18.5) and well separated (> 1'^) images, and to be accessible during the 2008B semester. The 
targets were also selected to have source redshifts < 2.8 so that L'-band observations correspond 
to rest-frame emission at > 1.0 /im. Information about the observations is given in Table [Tj 

All of our targets have previous ground- and space-based imaging. Results from single-epoch 
Hubble Space Telescope (HST) /, and H band photometry have been used in the past to derive 



lens models (e.g., Lehar et al. 2000). Among the various ground-based studies, optical monitoring 
programs are particularly useful to us because they help quantify the intrinsic variability in the 



source and the amplitude of variability induced by microlensing ( |Kochanek et al.||2006 ; Koptelova 



^Photometric monitoring can be used to quantify the variations, but has only been done for certain lenses ^ 


Eigen- 


brod et al.|2006||Kochanek et aL|2006| |Vuissoz et aL|2007 


|2008| iKoptelova et aL|2008| "Goicoechea et aL|2008 


). Four 


of our targets have been monitored at optical wavelengt 


hs: Q0142-100, HE 0435-1223, SBS 0909+523, and HE 


2149-2745. 


^ ht t p : / / www . cf a . harvard . edu / castles / 



et al.||2008| |Goicoechea et al.||2008D 



At wavelengths > 2 fim, thermal emission from Earth's atmosphere becomes an increasingly 
significant component of the flux recorded at the telescope. In order to remove the (often dominant) 
emission from the foreground sky, observations in the infrared are nodded between positions on and 
off source. Since our targets are much smaller than the field of view (image separation AO < 
compared with a 22^' field of view in L'), we decided to "nod on chip" by using a small 2x2, 
dither pattern. For each dither position, three other positions are available to construct sky model 
images. This strategy effectively doubles our observing efficiency by eliminating the need for off- 
source positions. In addition to providing sky subtraction images, dithering reduces the sensitivity 
of the observations to hot or bad pixels in the NIRI array. 

For K-band observations, single exposures of 20-40 s were selected to ensure that the bright 
quasar images (extrapolated from 1.6 to 2.2 /xm assuming fjy oc u' "Q-^; [Stern et al.||2005[ ) would not 



fill more than 50% of the detector well. In the V band, the sky brightness restricts exposure times 
to 1 s in order to prevent saturation. To accumulate the necessary integration time, the short V 
exposures were co-added in 30 s blocks at each dither position. 

Data reduction was conducted with standard IRAF routines provided by Gemini Observatory. 
The reduction proceeded as follows. First bad pixel masks and sigma clipping algorithms were 
applied to account for defective and hot pixels in the detector. The images were then flat fielded. 
For A^-band observations, sky brightness is subdominant and separate calibration flats were taken 
each night. Since the sky dominates L'-band observations, the images themselves are essentially sky 
flats and are used to construct flat field images. After flat fielding, the data were sky-subtracted 
using observations from other positions in the same dither pattern. sky images were limited 
to dither positions immediately before and after a given observation, since sky brightness can 



vary on the scale of minutes at 3.8 /xm (Glass 1999). Once each dither position was flat fielded 



sky subtracted, and sigma clipped, the final image was produced by combining the various dither 
positions. Figures [T] and [2] show the reduced K and U images for our six targets. 



3. Photometry 



We use the GALFIT image modeling package ( [Peng et aL]|2002 ) to extract the K and U fluxes 
of the lensed quasar images and (when detected) the lensing galaxy. GALFIT is valuable for general 
photometry, and specifically for separating the flux of the lensed images from that of the nearby 
lens galaxy, which is detected in our K-band images for all targets except SBS 0909+523. 



We follow a procedure similar to that of Sluse et al. (2008). We select the lensed quasar 



image farthest from the lens galaxy as our initial model for the point spread function (PSF)J^We 



^If the quasar images lie at similar distances from the lens galaxy, we select the brightest image. 
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then perform a simultaneous fit to all the lensed quasar images (modeled as scaled PSFs) and any 
galaxies present (modeled as PSF-convolved Sersic profiles). With this fit in hand, we subtract the 
galaxy (or galaxies) from the original image and then re-exact a new PSF model that contains less 
contamination from the lens galaxy. Using this cleaner PSF, we perform a new simultaneous fit. We 
compute the reduced Xr value for the residuals using regions defined by the PSF size for the quasar 
images, and regions that are 1, 2, or 3 times the PSF size for the lens galaxy. U Xr — ^ within these 
regions, we conclude that the procedure has converged. If not, we perform additional iterations, 
improving the model at each step, until convergence is achieved. Our A^-band images typically 
require 3 or 4 iterations, while our V images require only 1 or 2 since the lens galaxy is rarely 
detected in L^ GALFIT provides photometric errorbars, although each quasar image that serves 
as the PSF model is fit perfectly. For these images, we determine the uncertainty by computing 
the range in flux where Ax^ = 1 in the PSF region. (When the lensed quasar images are clearly 
distinct, we verify the derived flux values and uncertainties using aperture photometry.) Figure [s] 
shows an example of the GALFIT modeling procedure, and Table [2] gives the derived flux ratios. 

We calibrate the flux measurements using observations of standard stars taken immediately 
before or after science observations. We consider both aperture photometry and GALFIT modeling 
with the star itself providing the PSF model. Table |3] gives the calibrated photometry, along with 
a list of the standard stars we use. We report magnitudes in the Vega system. 
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Table 1. Targets and Observational Information 









Total Integration (min) 


Target ^lens ^source 


(") 


Date(s)t 


K 


V 


Q0142-100 0.491 2.719 


2.24 


8.21.2008 


3 


32 


SDSS 0246-0825 0.723 1.680 


1.04 


8.18.2008, 8.20.2008 4.8 


76 


HE 0435-1223 0.455 1.689 


2.42 


8.31.2008,12.22.2008 4.8 


64 


SDSS 0806+2006 0.573 1.540 


1.40 


12.22.2008 


2.4 


20 


SBS 0909+523 0.830 1.376 


1.17 


12.22.2008 


1.6 


2 


HE 2149-2745 0.603 2.033 


1.70 


12.22.2008 


3 


32 


'''Hawaii Standard Time 










Table 2. iiT and V Flux Ratios 




Target 


Image 


K 


L' 




Q0142-100 


A 


= 1.0 


= 1.0 






B 


0.128 + 0.002 


0.132 + 0.006 




SDSS 0246-0825 


A 


= 1.0 


= 1.0 






B 


0.258 + 0.015 


0.331 + 0.016 




HE 0435-1223 


A 


1.837 + 0.130 


1.706 + 0.129 






B 


1.271 + 0.073 


0.991 + 0.074 






C 


= 1.0 


= 1.0 






D 


0.851 + 0.049 


0.809 + 0.102 




SDSS 0806+2006 


A 


= 1.0 


= 1.0 






B 


0.406 + 0.030 


< 0.164 (3(t) 




SBS 0909+523 


A 


= 1.0 


= 1.0 






B 


0.973 + 0.028 


0.973 + 0.033 




HE 2149-2745 


A 


= 1.0 


= 1.0 






B 


0.280 + 0.006 


0.240 + 0.009 
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Table 3. K and L' Photometry 



Target 


Image 


K 


V 


Standard 


Q0142-100 


A 


14.61 + 0.01 + 0.009 


13.57 + 0.01 + 0.020 


1,2 - SA 92-342 {K,L') 




B 


16.84 + 0.02 + 0.009 


15.77 + 0.05 + 0.020 




SDSS 0246-0825 


A 


15.36 + 0.01 + 0.008 


13.71 + 0.01 + 0.008 


1 - GSPC S618-D {K) 




B 


16.83 + 0.06 + 0.008 


14.91 + 0.05 + 0.008 


2 - HD 22686 (L') 


HE 0435-1223 


A 


16.70 + 0.01 + 0.010 


15.00 + 0.02 + 0.010 


1,2 - HD 289907 {K,L') 




B 


17.10 + 0.02 + 0.010 


15.59 + 0.05 + 0.010 






C 


17.36 + 0.05 + 0.010 


15.58 + 0.05 + 0.010 






D 


17.68 + 0.05 + 0.010 


15.81 + 0.11 + 0.010 




SDSS 0806+2006 


A 


16.95 + 0.01 + 0.021 


15.25 + 0.01 + 0.020 


1,2 - GSPC P486-R {K,L') 




B 


17.93 + 0.08 + 0.021 


> 17.21 (3(t) 




SBS 0909+523 


A 


14.59 + 0.01 + 0.011 


13.30 + 0.01 + 0.013 


3 - HST P091-D [K) 




B 


14.57 + 0.03 + 0.011 


13.32 + 0.03 + 0.013 


2 - HD 84800 (L') 


HE 2149-2745 


A 


14.08 + 0.01 + 0.016 


14.23 ± 0.01 ± 0.006 


1 - EG 141 {K) 




B 


15.46 + 0.02 + 0.016 


15.78 + 0.04 + 0.006 


2 - GL 811.1 (L') 



Note. — Vega magnitudes of lensed images. The quoted uncertainties include contributions from the 
lensed images themselves (first) and the photometric standard (second). Col. 5 lists the standard star 
used for each passband. References: 1- Leggett et al. (2006), 2- |Leggett et al. (2003), 3- Persson et al. 

(pM). 
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Table 4. Previous Flux Ratio Measurements 



Target 


Filter 


B/A 


Ref. 


Q0142-100 


HST F555W 


0.127 + 0.009 


1 




HST F675W 


0.121 + 0.004 


1 




Johnson 


0.146 + 0.009 


2 




HST F814W 


0.146 + 0.016 


3 




HST F160W 


0.121 + 0.004 


1 


SDSS 0246-0825 


Sloan u 


0.319 + 0.033 


4 




Sloan g 


0.310 + 0.010 


4 




HST F555W 


0.483 + 0.010 


4 




Sloan r 


0.328 + 0.011 


4 




Sloan i 


0.340 + 0.016 


4 




HST F814W 


0.247 + 0.017 


4 




HST F160W 


0.291 + 0.010 


4 




Johnson H 


0.302 + 0.062 


4 




Johnson 


0.290 


4 


SDSS 0806+2006 


Johnson V 


0.581 + 0.012 


5 




Johnson R 


0.673 + 0.014 


5 




Johnson / 


0.581 + 0.012 


5 




0.39-1.10 /im 


0.7 


5 




0.45-0.87 /xm 


0.77 


6 




C III & Mg II BEL 


0.454 


6 




Johnson H 


0.649 + 0.013 


5 




Johnson H 


0.474 + 0.043 


6 




Johnson 


0.689 


5 


SBS 0909+523 


0.3-8 keV 


3.125 + 0.030 


7 




HST F555W 


0.444 + 0.077 


1 




Sloan 


0.551 + 0.012 


8 




HST F814W 


0.705 + 0.073 


1 




HST F160W 


0.887 + 0.029 


1 


HE 2149-2745 


Johnson B 


0.236 + 0.009 


9 




HST F555W 


0.209 + 0.015 


1 




Johnson V 


0.236 + 0.018 


9 




Johnson R 


0.231 + 0.006 


9 




HST F814W 


0.238 + 0.005 


1 
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Fig. 1.— K (left) and V (right) images of Q0142-100, SDSS 0246-0825, and HE 0435-1223 (top 
to bottom). Lensed images are labeled as A-D and lens galaxies are labeled as G or Gl, following 
nomenclature in the literature. Lens galaxies are clearly visible in the K band images. With the 
exception of Gl in SDSS 0246—0825, lens galaxies are not seen the V images. 




Fig. 2.— Similar to Figure ^ but for the lenses SDSS 0806+2006, SBS 0909+523, and HE 
2149-2745 (top to bottom). 
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Fig. 3. — A photometric model for the lens SDSS 0806+2006 produced by our GALFIT analysis of 
the K-band image in Figure [2j The top left panel shows the original data with the lens galaxy 
and quasar images labeled as G, A, and B, respectively. Image A provides the model PSF. The 
top right panel shows the model for the various components, while the bottom panel shows the 
residuals after model subtraction, which follow a Gaussian distribution and have reduced — 1. 
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4. Results 

In this section we use our new measurements to look for differences between the K and V flux 
ratios. In the next section we compare the K and V results with previous measurements at other 
wavelengths to discuss wavelength dependence more generally. 

• QOI42-IOO. The K and V flux ratios are Fb/Fa = 0.128 ± 0.002 and 0.132 ± 0.006, respec- 
tively, which are consistent within la. 

• SDSS 0246-0825. The flux ratios are Fb/Fa = 0.258 ± 0.015 and 0.331 ± 0.016 in the K 
and bands, respectively. The difference is inconsistent with measurement noise at > 99.9% 
confldence. 

• HE 0435— 1223. We consider flux ratios with respect to image C. For the A/C and D/C 
flux ratios, the K and V values are consistent within the errorbars. For the B/C flux ratio, 
however, there is a signiflcant difference with Fb/Fq = 1.271 ± 0.073 in K and 0.991 ± 0.074 
in L^ This discrepancy is signiflcant at > 99% confldence. 

• SDSS 0806+2006. In the K band the flux ratio is Fb/Fa = 0.406 ± 0.030. In the U band, 
however, we do not signiflcantly detect image B, with an upper limit of Fb/Fa < 0.164 (3a). 
This is the most striking wavelength dependence in our sample. 

• SBS 0909+523. We measure identical flux ratios of Fb/Fa = 0.973 in K and V. This nuh 
detection is interesting in comparison with known wavelength dependence in other bands (see 

• HE 2149-2745. The flux ratios are Fb/Fa = 0.280 ± 0.006 and Fb/Fa = 0.240 =b 0.009 in 
K and L', respectively, which differ at > 99.99% confldence. 



5. Discussion 

As discussed in ^ wavelength dependence in lens flux ratios (especially between V and shorter 
wavelengths) may signal effects of dark matter substructure, but it may also be associated with 
quasar variability, microlensing, or dust extinction. Here we compare our NIR flux ratio mea- 
surements (Table [2]) with previous measurements at other wavelengths (Table |4]) in an attempt 
to disentangle the various effects that can alter flux ratios across wavelengths. Even if we cannot 
uniquely identify the cause of wavelength dependence, our analysis highlights systems that have 
interesting features and warrant further study. To ease comparison, we plot in Figure [4] the flux 
ratios from Tables [2] and [4] between 0.4 — 3.8 /xm. 

QOI42— 100. Our K and flux ratios are compatible with each other and with single-epoch 



HST measurements in F555W, F675W, and F160W (Lehar et al. 2000). This is not surprising 



Table 4 — Continued 



Target 


Filter 




B/A 




Ref. 




HST F160W 




0.238 ± 0.009 




1 


Target 


Filter 


A/C 


B/C 


D/C 


Ref. 


HE 0435-1223 


HST F555W 


1.84 ±0.11 


1.08 ±0.09 


0.95 ± 0.06 


10 




KPNO i?t 


1.75 ±0.10 


1.00 ±0.04 


0.85 ± 0.05 


10 




Johnson 


1.69 ±0.06 


0.95 ±0.01 


0.82 ± 0.03 


11 




HST F814W 


1.69 ±0.04 


1.01 ±0.04 


0.82 ± 0.03 


10 




HST F160W 


1.57 ±0.05 


1.00 ±0.03 


0.79 ±0.03 


10 



^^Reported flux ratio is based on ground-based monitoring. Uncertainty includes 
scatter due to intrinsic and/or microlensing variability. 

•^Reported flux ratio is based on ground-based monitoring. Flux ratios are 
corrected for microlensing and intrinsic variabilty. 



Note. 

for our targets. References: 



This is a non-comprehensive list of existing flux ratio measurements 
Lehar et al.| ( |2000D , 2- |Koptelova et al.| ( [20Tol ), 3- 
CASTLES dhttp: / /www.cfa.harvard.edu/castles/D , 4-|Inada et al.| (|2005|), 5 



In- 



ada et al.| ( [20061 ), 6- |Sluse et aT1 ( |2008D , 7- |Dai fc Kochanek] ( [20091 ) , 8- [Goicoechea 



et al.[ ( [2008[ ), 9- |Lopez et aL| ( [1998[ ), 10- |Kochanek et al.| ( [2006D , ll- [Courbin et al. 
([20T0I). 
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because for a source redshift of Zg = 2.719 all bands up to L' probe rest-frame wavelengths < 1.0 /xm 



where the accretion disk is the main emission source (e.g., Honig et al. 2008). What is striking is 



that i?-band monitoring by Koptelova et al. (2010) yielded a mean flux ratio of Fb/Fa — 0.146 



and a scatter of 0.009 over 3-month observing windows in 2003-2005. That value, which agrees 
with the single-epoch F814W measurement from CASTLES, suggests that there is more intrin- 
sic/microlensing variability on several-year timescales than was apparent during the monitoring 
campaign. 

SDSS 0246— 0825. Previous observations in various filters and epochs point to significant 
variations with both time and wavelength. For instance, SDSS i and HST F814W measurements 



taken 20 months apart yielded Fb/Fa = 0.340 ± 0.016 and 0.247 ± 0.017, respectively (Inada 



et aL| 2005t CASTLES). At the same two epochs, SDSS g and HST F555W measurements gave 



Fb/Fa = 0.310 ± 0.010 and 0.483 ± 0.010. In the NIR, we find Fb/Fa = 0.258 ± 0.015 and 
0.331 ± 0.016 in K and L^ It is unclear whether the NIR discrepancy is produced by the same 
variations that affect the shorter wavelengths, or whether millilensing by dark matter substructure 
might be at play. To distinguish the possibilities, it would be interesting to use future multi- 
wavelength imaging to see whether optical and NIR flux ratios vary in similar ways. 

HE 0435—1223. The measured A/C and D/C flux ratios are consistent with being indepen- 
dent of wavelength over 0.5-3.8 /xm. Our measurement of the B/C flux ratio in V is consistent 



with measurements at other wavelengths (including the i?-band monitoring by Kochanek et al. 



2006[ ). However, the i^-band B/C flux ratio differs from the rest. This is curious because i^T-band 
corresponds to rest-frame optical emission (0.82 fim at Zg = 1.689), which should originate from 
the accretion disk. Using the standard thin-disk model of iShakura &; Sunyaev ( 1973[ ), we estimate 



the K source size to be only 4.7 larger than the R source size; for most plausible substructure 
models, such a difference is too small to account for the K-band anomaly (see iDobler &; Keeton 



2006). Combining that with the agreement between the R and flux ratios, we conjecture that 
the K-band observations in 2008 may have been affected by microlensing. Under this hypothesis, 
the L'-band observations were less sensitive to microlensing because the source was larger, while 
the earlier i?-band monitoring happened to catch a period of smaller microlensing variations. We 
study this interesting system in more detail in a separate paper (Fadely & Keeton, in prep.). 

SDSS 0806+2006. It is striking that we do not see image B in the L'-band observations 
{Fb/Fa < 0.164 at 3a), especially considering that Fb/Fa > 0.4 at shorter wavelengths. Changes 
in the flux ratio with both wavelength and time (see Table El) suggest that quasar variability and/or 



microlensing are at work in this lens (but not dust extinction; |Sluse et al. 2008). It is not clear 



though, whether either effect could explain the non-detection in L^ For Zg = 1.54 the U emission 
corresponds to rest-frame 1.5 /xm, so it presumably contains some contribution from the dusty 
torus, and thus ought to be less susceptible to intrinsic variations or microlensing than the shorter- 
wavelength bands. In any case, the wavelength dependence of the flux ratios is very interesting and 
should be investigated further with both modeling and follow-up observations. 
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SBS 0909-h523. There is significant differential extinction in this lens (iMotta et al.||2002|); 



our observations should further constrain dust models since they represent the longest wavelengths 



yet measured in this system. In fact, extrapolating the extinction models of Motta et al. (which 
were derived from shorter- wavelength measurements) yields predicted flux ratios of Fb/Fa = 0.909 
in K and 0.934 in V ^ whereas we measure 0.97 ± 0.03 in both bands. In terms of dark matter 
substructure, the agreement between our K and V flux ratios is interesting as a null detection. If 
the K (source rest-frame 0.9 /xm) and V (source rest-frame 1.6 /xm) sources have different effective 
sizes, there must be no substructure near the images with the right Einstein radius to produce a 
chromatic effect. Note that this does not rule out substructure in the lens — only substructure 
with the right mass at the right position. 

HE 2149—2745. The measured flux ratios are consistent with being independent of wavelength 
across 0.4-3.8 fim (Fb/Fa ^ 0.24), except that our K-band measurement (Fb/Fa = 0.280 ± 
0.006) formally disagrees with the others. However, the flux ratio is known to vary with recorded 
fluctuations at the level of 0.03 (Burud et al. 2002|). Since K band corresponds to rest-frame 
0.7 /xm, we conclude that the high i^-band flux ratio is likely due to intrinsic and/or microlensing 
variability. 



6. Conclusions 

We have used a multi- wavelength analysis of lens flux ratios to search for effects of dark matter 
substructure, along with quasar variability, microlensing, and dust extinction. We have presented 
new K and images of six lensed quasars, which were selected to have source redshifts < 2.8 so 
that observations probe rest-frame wavelengths > 1.0 /xm. Some of the flux should therefore 
originate from the extended torus of gas surrounding the central accretion disk, possibly providing 



the conditions for chromatic millilensing (see Dobler &; Keeton 2006) 



We find strong differences between the K and flux ratios for two lenses. In HE 0435—1223 
the U measurement of the B/C flux ratio is consistent with results at shorter wavelengths, but 
the K-band value is enhanced by ^30%. We argue in a separate analysis (Fadely & Keeton, 
in prep.) that the enhancement may be attributable to microlensing. To test that hypothesis, 
future observations are needed to look for variability in the X-band and compare it with variability 
at shorter wavelengths. In SDSS 0806+2006 we do not detect image B in V observations, with 
an upper limit of Fb/Fa < 0.164 (3cr), even though the flux ratio is Fb/Fa > 0.4 at shorter 
wavelengths. It would be very interesting to use deeper observations at U to find image B, to 
measure the flux ratio at other long wavelengths, and to combine all the observations with lens 
modeling to ascertain whether the U anomaly may be caused by dark matter substructure or some 
other interesting effect. 

We find a -30% difference between the K and V flux ratios for SDSS 0246-0825, and a 
smaller but statistically significant difference for HE 2149—2745. Multi-epoch observations of both 
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Fig. 4. — Flux ratios for our lens sample across 0.4 — 3.8 /im. New NIR (Table [2]) and existing 
measurements (Table [4]) are plotted with open and filled symbols, respectively. For HE 0435-1223 
we plot the flux ratios Fa/Fq^ Fb/Fq ^ 0.2, and Fd/Fq from top to bottom, with a thin dotted 
line added for clarity. 



systems have shown variations in the flux ratios at shorter wavelengths, which are presumably 
associated with quasar variability and/or microlensing. It is possible that the same phenomena 
affect the i^-band flux ratios and explain the NIR discrepancies. In that case, the flux ratios in 
V or other longer-wavelength bands would be useful for limiting effects of quasar variability and 
microlensing in order to probe dark matter substructure. Alternatively, it remains possible that 
chromatic millilensing causes the NIR discrepancies. To test both hypotheses, future observations 
are needed to quantify variability in multiple optical and NIR bands. 

Finally, we detect no difference between the K and L' flux ratios for Q0142— 100 and SBS 
0909+523. These results are interesting null detections, indicating that there is no substructure 
close to the images with masses in the right range to produce chromatic effects. Such results offer 
constraints on dark matter substructure that may be mild but nevertheless interesting for future 
statistical studies of substructure populations. 



We thank Gary Bernstein for helpful comments. This work was supported by NSF grant 



-17- 



AST-0747311. 

REFERENCES 

Agol, E., Gogarten, S. M., Gorjian, V., & Kimball, A. 2009, ApJ, 697, 1010 

Burud, L, Courbin, F., Magain, P., Lidman, C., Hutsemekers, D., Kneib, J., Hjorth, J., Brewer, 
J., Pompei, E., Germany, L., Pritchard, J., Jaunsen, A. O., Letawe, G., & Meylan, G. 2002, 
A&A, 383, 71 

Chartas, G., Kochanek, C. S., Dai, X., Poindexter, S., & Garmire, G. 2009, ApJ, 693, 174 
Chiba, M. 2002, ApJ, 565, 17 

Chiba, M., Minezaki, T., Kashikawa, N., Kataza, H., & Inoue, K. T. 2005, ApJ, 627, 53 

Courbin, F., Chantry, V., Revaz, Y., Sluse, D., Faure, C, Tewes, M., Eulaers, E., Koleva, M., 
Asfandiyarov, I., Dye, S., Magain, P., van Winckel, H., Coles, J., Saha, P., Ibrahimov, M., 
& Meylan, G. 2010, arXiv:1009.1473| 

Dai, X., & Kochanek, C. S. 2009, ApJ, 692, 677 

Dalai, N., & Kochanek, C. S. 2002, ApJ, 572, 25 

Diemand, J., Kuhlen, M., Madau, P., Zemp, M., Moore, B., Potter, D., & Stadel, J. 2008, Nature, 
454, 735 

Dobler, G., & Keeton, C. R. 2006, MNRAS, 365, 1243 

Eigenbrod, A., Courbin, F., Dye, S., Meylan, G., Sluse, D., Vuissoz, C, & Magain, P. 2006, A&A, 
451, 747 

Eliasdottir, A., Hjorth, J., Toft, S., Burud, L, & Paraficz, D. 2006, ApJS, 166, 443 

Glass, 1. S. 1999, Handbook of Infrared Astronomy, First Edition, 1st edn. (Cambridge University 
Press) 

Goicoechea, L. J., Shalyapin, V. N., Koptelova, E., Gil-Merino, R., Zheleznyak, A. P., & Ullan, A. 
2008, New Astronomy, 13, 182 

Honig, S. F., Prieto, M. A., & Beckert, T. 2008, A&A, 485, 33 

Inada, N., Buries, S., Gregg, M. D., Becker, R. H., Schechter, P. L., Eisenstein, D. J., Oguri, M., 
Castander, F. J., Hall, P. B., Johnston, D. E., Pindor, B., Richards, G. T., Schneider, D. P., 
White, R. L., Brinkmann, J., Szalay, A. S., & York, D. G. 2005, AJ, 130, 1967 



-18- 



Inada, N., Oguri, M., Becker, R. H., White, R. L., Gregg, M. D., Schechter, P. L., Kawano, Y., 
Kochanek, C. S., Richards, G. T., Schneider, D. P., Barentine, J. C, Brewington, H. J., 
Brinkmann, J., Harvanek, M., Kleinman, S. J., Krzesinski, J., Long, D., Neilsen, Jr., E. H., 
Nitta, A., Snedden, S. A., k York, D. G. 2006, AJ, 131, 1934 

Keeton, C. R., Gaudi, B. S., & Betters, A. O. 2003, ApJ, 598, 138 

Keeton, C. R., Gaudi, B. S., & Fetters, A. O. 2005, ApJ, 635, 35 

Kochanek, C. S., Morgan, N. D., Falco, E. E., McLeod, B. A., Winn, J. N., Dembicky, J., & 
Ketzeback, B. 2006, ApJ, 640, 47 

Koptelova, E. A., Artamonov, B. R, Bruevich, V. V., Burkhonov, O. A., & Sergeev, A. V. 2008, 
Astronomy Reports, 52, 270 

Koptelova, E., Oknyanskij, V. L., Artamonov, B. P., k Burkhonov, O. 2010, MNRAS, 401, 2805 

Leggett, S. K., Hawarden, T. G., Currie, M. J., Adamson, A. J. Carroll, T., Kerr, T. H., Kuhn, 
O. P., Seigar, M. S., Varricatt, W. P., feWold, T., 2003, MNRAS, 345, 144 

Leggett, S. K., Currie, M. J., Varricatt, W. P., Hawarden, T. G., Adamson, A. J., Buckle, J., 
Carroll, T., Davies, J. K., Davis, C. J., Kerr, T. H., Kuhn, O. P., Seigar, M. S., &Wold, T., 
2006, MNRAS, 373, 781 

Lehar, J., Falco, E. E., Kochanek, C. S., McLeod, B. A., Munoz, J. A., Impey, C. D., Rix, H., 
Keeton, C. R., & Peng, C. Y. 2000, ApJ, 536, 584 

Lopez, S., Wucknitz, O., & Wisotzki, L. 1998, A&A, 339, L13 

Maccio, A. V., Kang, X., Fontanot, F., Somerville, R. S., Koposov, S., & Monaco, P. 2010, MNRAS, 
402, 1995 

MacLeod, C. L., Kochanek, C. S., & Agol, E. 2009, ApJ, 699, 1578 
Madau, P., Diemand, J., & Kuhlen, M. 2008, ApJ, 679, 1260 
Metcalf, R. B., & Madau, P. 2001, ApJ, 563, 9 
Metcalf, R. B., & Zhao, H.-S. 2002, ApJ, 567, L5 

Minezaki, T., Chiba, M., Kashikawa, N., Inoue, K. T., & Kataza, H. 2009, ApJ, 697, 610 

Minezaki, T., Yoshii, Y., Kobayashi, Y., Enya, K., Suganuma, M., Tomita, H., Aoki, T., k Peterson, 
B. A. 2004, ApJ, 600, L35 

More, A., McKean, J. P., More, S., Porcas, R. W., Koopmans, L. V. E., k Garrett, M. A. 2009, 
MNRAS, 394, 174 



-19- 



Morgan, C. W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2010, ApJ, 712, 1129 

Motta, v., Mediavilla, E., Munoz, J. A., Falco, E., Kochanek, C. S., Arribas, S., Garcia-Lorenzo, 

B. , Oscoz, A., & Serra-Ricart, M. 2002, ApJ, 574, 719 

Moustakas, L. A., & Metcalf, R. B. 2003, MNRAS, 339, 607 

Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, Z., & Elitzur, M. 2008, ApJ, 685, 160 
Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H. 2002, AJ, 124, 266 

Persson, S. E., Murphy, D. C, Krzeminski, W., Roth, M., & Rieke, M. J., 1998, AJ, 116, 2475 

Pooley, D., Blackburne, J. A., Rappaport, S., t Schechter, P. L. 2007, ApJ, 661, 19 

Rowan-Robinson, M. 1995, MNRAS, 272, 737 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Sluse, D., Courbin, F., Eigenbrod, A., & Meylan, G. 2008, A&A, 492, L39 

Springel, V., Wang, J., Vogelsberger, M., Ludlow, A., Jenkins, A., Helmi, A., Navarro, J. F., Frenk, 

C. S., & White, S. D. M. 2008, MNRAS, 391, 1685 

Stern, D., Eisenhardt, P., Gorjian, V., Kochanek, C. S., Caldwell, N., Eisenstein, D., Brodwin, M., 
Brown, M. J. I., Cool, R., Dey, A., Green, P., Jannuzi, B. T., Murray, S. S., Pahre, M. A., 
& Willner, S. P. 2005, ApJ, 631, 163 

Vuissoz, C, Courbin, F., Sluse, D., Meylan, G., Chantry, V., Eulaers, E., Morgan, C, Eyler, M. E., 
Kochanek, C. S., Coles, J., Saha, P., Magain, P., & Falco, E. E. 2008, A&A, 488, 481 

Vuissoz, C, Courbin, F., Sluse, D., Meylan, G., Ibrahimov, M., Asfandiyarov, I., Stoops, E., 
Eigenbrod, A., Le Guillou, L., van Winckel, H., & Magain, P. 2007, A&A, 464, 845 

Wyithe, J. S. B., Webster, R. L., Turner, E. L., & Mortlock, D. J. 2000, MNRAS, 315, 62 



This preprint was prepared with the AAS IM^jX macros v5.2. 



